Genetic effects of heat stress on milk fatty acids in Brazilian Holstein cattle.

The present study aimed to estimate covariance components of milk fatty acids (FA) and to compare the genomic estimated breeding values under general and heat-stress effects. Data consisted of 38,762 test-day records from 6,344 Holstein cows obtained from May 2012 through January 2018 on 4 dairy herds from Brazil. Single-trait repeatability test-day models with random regressions as a function of temperature-humidity index values were used for genetic analyses. The models included contemporary groups, parity order (1-6), and days in milk classes as fixed effects, and general and thermotolerance additive genetic and permanent environmental as random effects. Notably, differences in heritability estimates between environments (general and heat stress) increased (0.03 to 0.06) for unsaturated FA traits, such as unsaturated, monounsaturated, and polyunsaturated, at higher heat-stress levels. In contrast, heritability estimated between environments for saturated FA traits, including saturated FA, palmitic acid (C16:0), and stearic acid (C18:0) did not observe significant differences between environments. In addition, our study revealed negative genetic correlations between general and heat-stress additive genetic effects (antagonistic effect) for the saturated FA, C16:0, C18:0, and C18:1, which ranged from -0.007 to -0.32. Spearman's ranking correlation between genomic estimated breeding values ranged from -0.27 to 0.99. Results indicated a moderate to strong interaction of genotype by the environment for most FA traits comparing a heat-stress environment with thermoneutral conditions. Our findings point out novel opportunities to explore the use of FA milk profile and heat-stress models.

Genome-wide association study for milk production traits in a Brazilian Holstein population.

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.

Human behavioral influences and milk quality control programs.

Mastitis is a major disease affecting the herds of dairy farmers worldwide. One of the indicators directly related to the widespread infection of this disease in herds is the bulk tank somatic cell count (BTSCC). Recent studies have shown that one of the risk factors associated with mastitis is the human factor. Therefore, understanding the influence of humans is essential to control and prevent the disease. The main goal of this study was to determine whether the motivations and barriers perceived by farmers could explain the variation in the BTSCC. This study was conducted at 75 dairy farms in southern Brazil. In the interviews with farmers, a survey based on Likert scale items was used to collect data. Structural equation models were used to explain the subjectivity in the ratio of observed variables and latent variables elucidating the possible causal relationships between the variables. The model indicated that some of the variation in the BTSCC can be explained by the farmer's behavior, which is elucidated by his/her motivations and barriers. The correlations between motivations and the BTSCC and between barriers and the BTSCC were positive. These findings suggest that variations in the BTSCC can be explained by the motivations and barriers perceived by farmers and that the Fogg Behavior Model used in this study can be used to explain how human behaviors influence mastitis control. This study also indicates that consulting companies focused on improving milk quality should pay attention to the human factor to reduce these barriers.

Do different standard plate counting (IDF/ISSO or AOAC) methods interfere in the conversion of individual bacteria counts to colony forming units in raw milk?

AIMS: This study aimed to establish the correlation between individual bacterial count (IBC) obtained by flow cytometry and the number of colony forming units (CFU) determined by standard plate count (SPC) in raw milk using two different reference methodologies: the methodology of the International Dairy Federation (IDF) - International Organization for Standardization (ISO) 4833, incubation for 72 h at 30°C and the methodology of the Association of Official Analytical Chemists (AOAC), incubation for 48 h at 35°C. METHODS AND RESULTS: For this, 100 bovine milk samples (80 ml) from different farms were collected in a sterile bottle and maintained refrigerated at 4°C and were delivered to the laboratory. In the laboratory, the samples were divided into two vials of 40 ml each. Then, half of the vials were forwarded for the SPC analysis, and the other half were analysed using the equipment BactoScan FC. The analyses by flow cytometry and SPC were performed at the same time (maximum deviation of +/- 1 h). To transform the data from IBC ml(-1) to CFU ml(-1) (IDF or AOAC methodology), a standard linear regression equation was used, as recommended by IDF/ISO-196. The difference between the reference methodologies affects the equation that transforms IBC into CFU and therefore the accuracy of the results. The results estimated by the equation using the ISO 4833 methodology were on average 0·18 log units higher than the results estimated using the equation using the AOAC methodology. CONCLUSIONS: After the comparison of the methodologies, it was concluded that there is an impact of the reference methodologies on the conversion of the results from IBC to CFU. SIGNIFICANCE AND IMPACT OF THE STUDY: Depending on the methodology adopted by each laboratory or country, there may not be equivalence in the results. Hence, the laboratories specialized in milk quality analysis that have changed their methodology for analysis, passing from the MAPA (AOAC) methodology to the IDF standard, need to develop new conversion equations to make their results in CFU ml(-1) consistent regardless of the reference methodology.

Genetic parameters for milk fatty acids, milk yield and quality traits of a Holstein cattle population reared under tropical conditions.

Information about genetic parameters is essential for selection decisions and genetic evaluation. These estimates are population specific; however, there are few studies with dairy cattle populations reared under tropical and sub-tropical conditions. Thus, the aim was to obtain estimates of heritability and genetic correlations for milk yield and quality traits using pedigree and genomic information from a Holstein population maintained in a tropical environment. Phenotypic records (n = 36 457) of 4203 cows as well as the genotypes for 57 368 single nucleotide polymorphisms from 755 of these cows were used. Covariance components were estimated using the restricted maximum likelihood method under a mixed animal model, considering a pedigree-based relationship matrix or a combined pedigree-genomic matrix. High heritabilities (around 0.30) were estimated for lactose and protein content in milk whereas moderate values (between 0.19 and 0.26) were obtained for percentages of fat, saturated fatty acids and palmitic acid in milk. Genetic correlations ranging from -0.38 to -0.13 were determined between milk yield and composition traits. The smaller estimates compared to other similar studies can be due to poor environmental conditions, which may reduce genetic variability. These results highlight the importance in using genetic parameters estimated in the population under evaluation for selection decisions.

Efeito da contagem de células somáticas sobre o rendimento e a composição físico-química do queijo muçarela / The effect of somatic cell count on yield and physico-chemical composition of Mozzarella cheese

Determinou-se o efeito da contagem de células somáticas do leite sobre o rendimento e a composição físico-química do queijo muçarela. Foram selecionadas vacas com contagem de células somáticas (CCS) <200.000 células/mL; CCS >200.000 a <400.000 células/mL; CCS >400.000 células/mL a <750.000 células/mL e CCS >750.000 células/mL e que não receberam tratamento com antimicrobianos nos dias que antecederam a obtenção do leite utilizado no processamento. Os queijos foram avaliados por meio das análises de pH, acidez, umidade, proteína, gordura, extrato seco total e desengordurado. Os queijos elaborados do leite com alta CCS >750.000 cél./mL apresentaram menor teor de proteína, maior umidade e menor rendimento industrial. Conclui-se que a elevação da CCS cursou com alterações na qualidade do queijo muçarela, o que reflete em sua composição nutricional e microbiológica, pois ocorre um incremento na atividade de água, fator intrínseco indispensável para o crescimento microbiano e consequente redução da vida de prateleira...

The effect of somatic cell count (SCC) on the yield and composition of mozzarella cheese was evaluated. Cows with different levels of SCC in their milk (<200.000 cells/mL; >200.000 to <400.000 cells/mL; >400.000 cells/mL to <750.000 cells/mL and >750.000 cells/mL) were selected. The animals had no antimicrobial treatment prior to or on milk sampling day. The mozzarella cheese quality was evaluated analyzing pH, acidity, moisture and protein content, fat, total solids and nonfat solids. The trial was replicated four times and the experimental design was analyzed through randomized blocks. The mozzarella cheese from milk with high SCC (>750.000 cells/mL) had lower protein content, higher moisture content and lower industrial yield, compared to cheese from milk with lower levels of SCC. High numbers of SCC impaired mozzarella cheese quality, which was reflected in its nutritional and microbiological composition. We observed an increase in water activity, an intrinsic and essential factor for microbial growth...

Short communication: Principal components and factor analytic models for test-day milk yield in Brazilian Holstein cattle.

A total of 46,089 individual monthly test-day (TD) milk yields (10 test-days), from 7,331 complete first lactations of Holstein cattle were analyzed. A standard multivariate analysis (MV), reduced rank analyses fitting the first 2, 3, and 4 genetic principal components (PC2, PC3, PC4), and analyses that fitted a factor analytic structure considering 2, 3, and 4 factors (FAS2, FAS3, FAS4), were carried out. The models included the random animal genetic effect and fixed effects of the contemporary groups (herd-year-month of test-day), age of cow (linear and quadratic effects), and days in milk (linear effect). The residual covariance matrix was assumed to have full rank. Moreover, 2 random regression models were applied. Variance components were estimated by restricted maximum likelihood method. The heritability estimates ranged from 0.11 to 0.24. The genetic correlation estimates between TD obtained with the PC2 model were higher than those obtained with the MV model, especially on adjacent test-days at the end of lactation close to unity. The results indicate that for the data considered in this study, only 2 principal components are required to summarize the bulk of genetic variation among the 10 traits.

Random regression models using different functions to model test-day milk yield of Brazilian Holstein cows.

We analyzed 152,145 test-day records from 7317 first lactations of Holstein cows recorded from 1995 to 2003. Our objective was to model variations in test-day milk yield during the first lactation of Holstein cows by random regression model (RRM), using various functions in order to obtain adequate and parsimonious models for the estimation of genetic parameters. Test-day milk yields were grouped into weekly classes of days in milk, ranging from 1 to 44 weeks. The contemporary groups were defined as herd-test-day. The analyses were performed using a single-trait RRM, including the direct additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. The mean trend of milk yield was modeled with a fourth-order orthogonal Legendre polynomial. The additive genetic and permanent environmental covariance functions were estimated by random regression on two parametric functions, Ali and Schaeffer and Wilmink, and on B-spline functions of days in milk. The covariance components and the genetic parameters were estimated by the restricted maximum likelihood method. Results from RRM parametric and B-spline functions were compared to RRM on Legendre polynomials and with a multi-trait analysis, using the same data set. Heritability estimates presented similar trends during mid-lactation (13 to 31 weeks) and between week 37 and the end of lactation, for all RRM. Heritabilities obtained by multi-trait analysis were of a lower magnitude than those estimated by RRM. The RRMs with a higher number of parameters were more useful to describe the genetic variation of test-day milk yield throughout the lactation. RRM using B-spline and Legendre polynomials as base functions appears to be the most adequate to describe the covariance structure of the data.

Sazonalidade do teor de proteína e outros componentes do leite e sua relação com programa de pagamento por qualidade / Seasonability of protein and others milk components related with quality payment program

Avaliou-se o efeito da variação da qualidade do leite mediante acompanhamento de amostras provenientes de tanques refrigeradores. As 2.970 amostras foram obtidas de outubro/2005 a setembro/2006. Foram realizadas análises de composição, contagem de células somáticas (CCS) e contagem bacteriana total. As condições climáticas foram avaliadas pela temperatura e umidade do ar coletadas em estações meteorológicas distribuídas em todo o estado de São Paulo. No outono, ocorreram os maiores valores para gordura e proteína do leite, 3,65 e 3,21 por cento, respectivamente. Para contagem de células somáticas, foi observado aumento no verão e contagem bacteriana na época da primavera. Foi identificada a influência da época do ano sobre o preço do leite pago dentro de um programa, com valor de bonificação máxima de 7,9 por cento para o outono. Um dos mais importantes itens relacionados com a penalização foi a CCS, a qual pode interferir na qualidade dos produtos finais por meio de alteração da proporção de caseína/proteínas do soro. O efeito da CCS e o da época do ano devem ser considerados na formulação das tabelas de programas de pagamento.

This study evaluated the variation of quality of bulk tank milk using 2,970 samples, collected from October 2005 to September 2006. Milk composition analyses, somatic cell count (SCC), and total bacterial count were performed. Climate conditions were evaluated by temperature and relative air humidity collected at meteorological stations throughout São Paulo state. Autumn season was related to the highest values for protein and fat, 3.65 and 3.21 percent, respectively. For somatic cell count, it was observed an increase in the count in the summer and for bacterial count in the spring. The influence of season on milk price was observed for dairy producers in a program, with maximum bonus of 7.9 percent in the autumn. One of the most important parameters related to penalty for lower quality of milk was SCC, since it can affect the quality of some dairies by the alteration of casein/whey protein proportion. Therefore, the effects of SCC and season should be considered into the process of formulation of tables for milk payment program.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models.

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Short communication: evaluation of an on-farm test to estimate somatic cell count.

The objective of this study was to compare the results of an on-farm test, named Somaticell, with results of electronic cell counting and for milk somatic cell count (SCC) among readers. The Somaticell test correctly determined the SCC in fresh quarter milk samples. Correlation between Somaticell and electronic enumeration of somatic cells was 0.92 and kappa coefficient 0.82. Using a threshold of 205,000 cells/mL, the sensitivity and specificity for determination of intramammary infections were 91.3 and 96.0%, respectively. The SCC was greater for milk samples from which major mastitis pathogens were recovered. Minor variation among readers was observed and most likely associated with the mixing procedure. However, the final analysis indicated that this variation was not significant and did not affect the amount of samples classified as having subclinical mastitis. The on-farm test evaluated in this study showed adequate capacity of determining SCC on quarter milk samples and may be considered as an alternative for on-farm detection of subclinical mastitis.

Genetic parameters for test day milk yields of first lactation Holstein cows by random regression models.

Covariance components for test day milk yield using 263 390 first lactation records of 32 448 Holstein cows were estimated using random regression animal models by restricted maximum likelihood. Three functions were used to adjust the lactation curve: the five-parameter logarithmic Ali and Schaeffer function (AS), the three-parameter exponential Wilmink function in its standard form (W) and in a modified form (W*), by reducing the range of covariate, and the combination of Legendre polynomial and W (LEG+W). Heterogeneous residual variance (RV) for different classes (4 and 29) of days in milk was considered in adjusting the functions. Estimates of RV were quite similar, rating from 4.15 to 5.29 kg2. Heritability estimates for AS (0.29 to 0.42), LEG+W (0.28 to 0.42) and W* (0.33 to 0.40) were similar, but heritability estimates used W (0.25 to 0.65) were highest than those estimated by the other functions, particularly at the end of lactation. Genetic correlations between milk yield on consecutive test days were close to unity, but decreased as the interval between test days increased. The AS function with homogeneous RV model had the best fit among those evaluated.

Análise econômica de esquemas alternativos de controle leiteiro / Economic analysis of alternative schemes of milk recording

Compararam-se os custos de seis esquemas de controle leiteiro (CL), estabelecidos com base em diferentes freqüências de controles (com possibilidade de controle trimestral supervisionado e aproveitamento do CL da fazenda), no número de ordenhas e na realização de análises qualitativas (composição do leite e/ou contagem de células somáticas). Computaram-se os percentuais que os custos totais do CL representariam na receita bruta mensal do leite (CL/RB) e na receita líquida mensal (CL/RL). O esquema tradicional (freqüência mensal) apresentou o mais elevado custo mensal e os maiores CL/RB e CL/RL. Esquemas com intervalo maior entre controles e amostragem em ordenhas alternadas resultaram na redução dos custos com diárias e quilometragem. O custo mensal do CL variou de 0,68% a 1,8% sobre a receita bruta do leite e de 6,6% a 17,0% sobre a receita líquida.

Indicadores de desempenho em rebanho da raça Holandesa: curvas de crescimento e altura, características reprodutivas, produtivas e parâmetros genéticos / Productive indicator in a Holstein dairy herd: growth and height curves, reproductive and productive traits and genetic parameters

Foram analisados dados de peso, altura e desempenho produtivo e reprodutivo de 1.272 novilhas da raça Holandesa, nascidas no período de 1990 a 1998. Estudaram-se o desenvolvimento em peso e altura do nascimento até o primeiro parto, as características peso ao primeiro parto (PP), idade ao primeiro parto (IP), produção de leite na primeira lactação (PL), primeiro intervalo de partos (IEP), e os parâmetros genéticos pertinentes. As equações para as curvas de peso e altura, estimadas pelo modelo proposto por von Bertalanffy, foram: Yi= 802,10 (1 - 0,630EXP (-0,0022t)) e Yi= 1,40 (1 - 0,208EXP (-0,0038t)), respectivamente. As novilhas atingiram o primeiro parto com altura média de 140cm e peso de 602kg aos 27 meses de idade. A produção total de leite na primeira lactação foi de 8.026kg e o primeiro IEP 420 dias. Na obtenção dos componentes de (co)variância utilizou-se o método da máxima verossimilhança restrita livre de derivada. As estimativas de herdabilidade foram 0,31, 0,35 e 0,34, respectivamente, para PP, IP e PL, e 0,10 para IEP, evidenciando variância genética aditiva média para as três primeiras características. As correlações genéticas apresentaram pequena magnitude para PP e PL, PL e IP, e IEP e PL. A correlação entre IEP e PP foi 0,49; entre IEP e IP o valor foi negativo, -0,69, indicando provável antagonismo entre as duas últimas características.

Prediçäo da taxa de gestaçäo de novilhas da raça Nelore acasaladas com um ano de idade / Prediction of pregnancy rate of Nellore heifers mated at yearling age

Os objetivos deste trabalho foram quantificar o relacionamento entre variáveis observadas antes da estaçäo de monta e a fertilidade de 321 novilhas e usar tais variáveis no desenvolvimento de equaçöes para predizer a fertilidade das novilhas usando-se metodologia de regressäo logística. As variáveis foram: pesos ao nascer, à desmama e a um ano de idade, ganhos de peso diários na pré e pós-desmama e no total, taxas de crescimento relativo(TCR) na pré e pós-desmama e no total, idade da novilha no início da estaçäo de monta (IDCOB) e idade da vaca ao parto. As novilhas que conceberam eram 18 dias mais velhas do que as que näo conceberam na estaçäo. As novilhas que ficaram gestantes eram 1,1kg mais leves ao nascer do que as que falharam em conceber. As que tiveram diagnóstico de gestaçäo positivo apresentaram, em média, 0,01(porcento) a mais nas taxas de crescimento relativo pré-desmama e total do que as que falharam em conceber. Novilhas que ficaram gestantes eram filhas de vacas mais jovens ao parto e o modelo completo explicou 26,6(porcento) da variaçäo da fertilidade delas. A combinaçäo dos efeitos lineares da IDCOB e da TCR total foi o melhor modelo para predizer a fertilidade das novilhas

Correlaçäo entre contagens de células somáticas determinadas no "Somacount 300" utilizando-se os detergentes RBS 35 e Extram MA 01 alcalino / Correlation between milk somatic cell counting of samples prepared with two detergents, and analyzed by the Somacount

The objective of this work was to compare the values of milk somatic cell counts (SCC) using the Somacount 300, and two detergents for sample preparation. the equipment was calibrated using the detergent RBS 35. Forty independent milk samples were analyzed. Each sample was randomly divided into two subsamples, and each subsample was analyzed after preparation with either RBS 35 or Extram MA 01 Alkaline detergents. the correlation coefficient between the two methods was high (0.998). This finding suggests the possibility of substitution of RBS 35 for an economical alternative

Effect of the combination of monensin and isoacids on rumen fermentation in vitro.

Effects of isoacids, monensin, or a combination of them on fermentation by mixed rumen bacteria were investigated using a continuous culture technique. The culture was allowed to stabilize for 4 d before treatments were imposed. Comparisons between treatments were made on d 11 and 12 of the culture. Isoacids (equal proportions of isobutyric, 2-M-butyric, isovaleric, and valeric acids) at 15 mg/dl of culture media increased acetate (6.17 vs. 5.48 meq/dl) and total VFA production (8.93 vs. 7.87 meq/dl) compared with that of controls. Monensin at 150 micrograms/dl reduced acetate (3.74 vs. 6.02 meq/dl) and VFA (6.84 vs. 8.54 meq/dl) but increased propionate (2.28 vs. 1.74 meq/dl) relative to control. The combination of isoacids and monensin increased acetate relative to monensin alone (5.24 vs. 3.74 meq/dl) but did not alter the effect of monensin on propionate concentration (2.32 vs. 2.28 meq/dl). It is concluded that monensin decreases acetate production by 35% and when isoacids are added to the cultures containing monensin, acetate production is restored.